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Abstract. The purpose of this work is to study the structure and evolution of white dwarf stars with helium-rich 
atmospheres (DB) in a self-consistent way with the predictions of time-dependent element diffusion. Specifically, 
we have considered white dwarf models with stellar masses in the range O.6O-O.85Af0 and helium envelopes 
with masses from 10 -2 to 10 _4 M». Our treatment of diffusion, appropriate for multicomponent gases, includes 
gravitational settling and chemical and thermal diffusion. OPAL radiative opacities for arbitrary metallicity 
and carbon-and oxygen-rich compositions are employed. Emphasis is placed on the evolution of the diffusion- 
modeled double-layered chemical structure. This structure, which is characterized by a pure helium envelope atop 
an intermediate remnant shell rich in helium, carbon and oxygen, is expected for pulsating DB white dwarfs, 
assuming that they are descendants of hydrogen-deficient PG1159 post-AGB stars. We find that, depending on 
the stellar mass, if DB white dwarf progenitors are formed with a helium content smaller than ~ 10 _3 M«, a 
single-layered configuration is expected to emerge during the DB pulsation instability strip. We also explore the 
consequences of diffusively evolving chemical stratifications on the adiabatic pulsational properties of our DB white 
dwarf models. In this context, we find that the evolving shape of the chemical profile translates into a distinct 
behaviour of the theoretical period distribution as compared with the case in which the shape of the profile is 
assumed to be fixed during the evolution across the instability strip. In particular, we find that the presence of a 
double-layered structure causes the period spacing diagrams to exhibit mode-trapping substructures. Finally, we 
extend the scope of the calculations to the domain of the helium-rich carbon-contaminated DQ white dwarfs. In 
particular, we speculate that DQ white dwarfs with low detected carbon abundances would not be descendants 
of the PG1159 stars. 
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1. Introduction 

White dwarf stars with helium-rich atmospheres (DB) 
constitute about the 20 % of all observed white dwarfs. 
The majority of these stars are thought to be the result of 
a very late thermal pulse experienced by post-asymptotic- 
giant-branch (post-AGB) progenitors during their early 
cooling phase (the born-again scenario; see Iben et al. 
1983). As a result of such a pulse, most of the residual hy- 
drogen envelope is completely burnt, and the star returns 
to the red giant region and then to the planetary nebula 
regime at high effective temperature (T c g ) as a helium- 
burning object. Evolutionary calculations which incorpo- 
rate diffusive overshooting (Herwig et al. 1999) show that 
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the occurrence of a very late thermal pulse leads to the for- 
mation of hydrogen-deficient post-AGB stars with surface 
abundances of helium, carbon and oxygen. This prediction 
is in good agreement with the surface abundance pattern 
observed in most PG1159 stars or the WR-type central 
stars of planetary nebulae (Drcizlcr & Hcber 1998), the 
hot and hydrogen-deficient remnants of post-AGB stars 
generally believed to be the immediate progenitors of most 
DB white dwarfs. 

The detection of pulsations in some DB white dwarfs 
(in the temperature range of ~ 27000 — 20000K) has al- 
lowed researchers to provide independent valuable con- 
straints on the internal structure and evolution of these 
stars, as well as on their progenitors. Indeed, asteroseis- 
mological techniques have become a powerful tool for in- 
ferring fundamental properties such the white dwarf mass, 
effective temperature, core composition and helium layer 
mass (see Bradley & Winget 1994; Metcalfe et al. 2000 and 
Metcalfe et al. 2001). In particular, Metcalfe et al. 2001 
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(see also Metcalfe et al. 2002) have recently applied DB 
white dwarf asteroseismology to place constraints on the 
12 C(a,7) 16 reaction rate from inferences for the abun- 
dance of central oxygen in the pulsating DB GD 358. It is 
worth noting that such studies are based on stellar models 
for pulsating DBs with a pure helium envelope atop the 
carbon-oxygen core (single-layered envelope). 

However, recent theoretical evidence makes past infer- 
ences about core composition of DBs and the 12 C(a, 7) 16 
reaction rate somewhat uncertain. Indeed, Fontaine & 
Brassard (2002) have shown that the diffusion-built-up, 
double-layered chemical structure expected for most of the 
pulsating DBs leads to a distinct theoretical period dis- 
tribution from that predicted by stellar models in which 
a single-layered configuration is assumed. The presence 
of a double-layered envelope in DB white dwarfs has 
been found in evolutionary calculations including time- 
dependent element diffusion by Dehner & Kawaler (1995) 
and Gautschy & Althaus (2002), assuming that DB white 
dwarfs descend from PG 1159 stars. Indeed, such calcu- 
lations show that by the time the domain of the variable 
DBs is reached, models are characterized by two different 
chemical transition zones (double-layered configuration). 
In fact, above the carbon-oxygen core, there exists an en- 
velope consisting of an intershell region rich in helium, 
carbon and oxygen, the relics of the short-lived mixing 
episode occurred during the last helium thermal pulse, 
and an overlying pure helium mantle resulting from grav- 
itational settling of carbon and oxygen. Finally, DB as- 
teroseismological fittings incorporating both the double- 
layered envelope feature expected from time-dependent 
diffusion calculations and adjustable carbon-oxygen cores 
have recently been presented by Metcalfe et al. (2003) 
for a wide range of helium contents and stellar masses. 
Despite their models yielding significantly better fits to 
observations, the derived stellar parameters for some fit- 
tings led them to conclude that double-layered models 
with adjustable carbon/oxygen cores may not be entirely 
appropriate to explain the observations. 

In view of these concerns, this paper is aimed at explor- 
ing the evolution of DB white dwarfs in a self-consistent 
way with the predictions of time-dependent element diffu- 
sion. Emphasis is placed on the chemistry variations along 
the evolutionary stages during which these stars are ex- 
pected to pulsate. Specifically, we hope to address the fol- 
lowing question: Is it possible that the diffusion-induced, 
double-layered structure could be altered by further chem- 
ical evolution to such an extent that a single-layered con- 
figuration emerges before the star reaches the red edge of 
the instability strip?. In view of the lack of previous cal- 
culations of the mass-dependence of time-dependent diffu- 
sion profiles, we judge that this question is an important 
one to be answered. Additionally, we extend the scope 
of the paper and explore the consequences of diffusively 
evolving chemical stratifications for the adiabatic pulsa- 
tional properties of the models. Finally, the implications 
of our results for the carbon surface abundance in some 
helium-rich DQ white dwarfs are explored. 



2. Computational details and input physics 

The results presented in this work have been obtained 
with the DB white dwarf evolutionary code developed at 
La Plata Observatory. Except for minor modifications, the 
code is basically that described in Gautschy & Althaus 
(2002) and references cited therein. In particular, micro- 
physics includes an updated version of the equation of 
state of Magni & Mazzitelli (1979), OPAL radiative opac- 
ities for arbitrary metallicity (Iglesias & Rogers 1996) in- 
cluding carbon- and oxygen-rich compositions, and up-to- 
date neutrino emission rates and conductive opacities (see 
Althaus et al. 2002). In particular, opacities for various 
mctallicities are required because of the metallicity gra- 
dient that develops in the envelopes as a result of gravi- 
tational settling and also in the outer layers of DQ white 
dwarfs because of dredge-up episodes (see next section). 
In this work, convection is treated in the framework of the 
mixing length theory as given by the ML2 parametrization 
(see Tassoul et al. 1990). 

The evolution of the chemical abundance distribution 
caused by diffusion processes has been fully accounted for 
by means of a time-dependent, finite-difference scheme 
that solves elemental continuity equations. To compute 
the white dwarf evolution in a self-consistent way with 
the predictions of element diffusion, such a scheme is cou- 
pled to the evolutionary code. In our treatment of time- 
dependent diffusion we have considered gravitational set- 
tling and chemical and thermal diffusion for the nuclear 
species 4 He, 12 C and 16 0. Radiative levitation, which is 
important in connection with surface composition of hot 
white dwarfs has been neglected. Specifically, we have 
adopted the treatment of Burgers (1969) for multicom- 
ponent gases. There, diffusion velocities Wi and residual 
heat flows satisfy the set of equations 
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Here, pi, pi, rii, Zi and mj means, respectively, the par- 
tial pressure, mass density, number density, mean charge 
and mass for species i (N means the number of ionic 
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Fig. 1. Abundance by mass of 4 He and 12 C (solid and 
dashed line, respectively) as a function of the outer mass 
fraction q for 0.60-M Q DB models at T cS = 27300^. Thick 
(thin) lines correspond to the situation in which thermal 
diffusion is considered (neglected). Clearly, the inclusion 
of thermal diffusion leads to thicker pure helium mantles. 

species plus electron). T, ks and VT are the temperature, 
Boltzmann constant and temperature gradient, respec- 
tively. The unknown variables are Wi and Ti (for ions and 
electrons). The electric field E has also to be determined. 



z\ - and z\'-) are from 



The resistance coefficients (K,, 

Paquette at al. (f986). For more details, see Gautschy & 
Althaus (2002). ft is worth mentioning that a recent study 
of the 3 He diffusion in DB white dwarfs has been presented 
by Montgomery et al. (200f ). 

The shape of the chemical profile in the envelope of 
a white dwarf is determined by the competition between 
basically partial pressure gradients, gravity and induced 
electric fields. In order to get an insight into the essential 
physics driving element diffusion and to understand some 
of the results to be presented later we will assume for 
a moment a simplified situation in which only two ionic 
species with mean charge Z\ and Z 2 and atomic mass 
number A\ and A 2 are present. In addition, we assume 
the electron to have zero mass and neglect the contribution 
from thermal diffusion. It is then straightforward to show 
from Eq. (I) that the diffusion velocity satisfies 

K t ( A * A A 
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where w\ is the diffusion velocity for specie I < 
means that the element sinks into the star) and / = 



(f + Aini/A 2 n 2 ){niZi +n 2 Z 2 ) /n 1 Z 1 n 2 Z 2 . In Eq. (3), m H 
and g are the hydrogen-atom mass and gravity, respec- 
tively. In deriving Eq. (3) we have assumed an ideal gas 
for ions and use the condition for no net mass flow rel- 
ative to the centre of mass. The first term on the right- 
hand side of Eq. (3) takes into account the contribution 
of the gravitational settling and the influence of the in- 
duced electric field, and the third and fourth terms refer 
to the chemical diffusion contribution. We apply Eq. (3) 
to a mixture of 4 He and 12 C. Note that the ionization 
state of elements strongly affects the diffusion velocity. In 
particular, in the non-degenerate outer layers, 12 C is not 
fully ionized (for instance, in our T c s = 27300 K, 0.6-M Q 
model, 12 C becomes completely ionized at a mass depth 
of 10 _8 M* below the stellar surface), and as a result 4 He 
diffuses upwards very rapidly driven mainly by gravity 
(first term in Eq. 3 is dominant). As cooling proceeds, 
this leads to a thickening pure helium mantle. In deeper, 
more degenerate layers, gravitational settling turns out to 
be less operative. Indeed, 12 C (and 4 He) becomes fully 
ionized and the gravity term vanishes. There, diffusion is 
driven essentially by chemical gradients, thus forcing he- 
lium to penetrate downwards to hotter layers. Note that 
this is quite different from the situation for a mixture of 
1 H and 4 He. In that case, the gravity term never vanishes 
and, at large degeneracy, it takes over chemical diffusion, 
thus halting the inward diffusion of *H and leading to a 
complete separation of the two species at advanced ages. 

It is important to note that our diffusion treatment 
considers thermal diffusion. This is by no means a negligi- 
ble contribution to the diffusion processes. Because ther- 
mal diffusion acts in the same direction as gravitational 
settling, the neglect of thermal diffusion underestimates 
the rate at which carbon and oxygen diffuse downwards. 
Hence, at a given T c g value, we expect a thicker pure he- 
lium mantle when thermal diffusion is taken into account. 
This expectation is indeed borne out by Fig. ^ which il- 
lustrates the helium and carbon abundance distribution 
in terms of the outer mass fraction q (q = 1 — M r /M„) for 
a 0.6-Mq DB white dwarf model at T cS = 27300 K, i.e. 
near the blue edge of the DB instability strip. Note that 
the mass embraced by the pure helium mantle amounts to 
2.5 x 10~ 6 M*, as compared with the value of 6 x f §~ 7 
when thermal diffusion is neglected. 

With regard to the initial stellar models needed to start 
our cooling sequences, they were obtained by means of 
the same artificial procedure as described in Gautschy & 
Althaus (2002). In particular, such starting models cor- 
respond to hydrogen-deficient post-AGB configurations 
with an uniform envelope chemical composition of helium, 
carbon and oxygen (with abundances of 0.42, 0.36 and 0.22 
by mass, respectively). As mentioned in the introduction, 
the envelope composition is the result of a short-lived mix- 
ing event occurred during the late thermal pulse experi- 
enced by a post-AGB remnant. The chemical composition 
of the core is that of Salaris et al. (1997). Specifically, 
we have considered DB white dwarf models with stellar 
masses of 0.60, 0.65, 0.70, 0.75 and 0.85 M . The mass 
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of the helium envelope (Mh c ) is constrained by the the- 
ory of post-AGB evolution to be in the range 10~ 2 to 10 -3 
M* , in good agreement with carbon dredge- up simulations 
in DQ white dwarfs, the cooler descendants of DBs, by 
MacDonald et al. (1998) for a typical white dwarf mass of 
0.6 Mq. In this work, we considered the following values 
for M He : 8 x 1CT 3 ,9 x 1CT 4 and 1CT 4 M*. 

In what follows, we describe the main results of our 
calculations. We want to mention that this paper is not 
aimed at performing asteroseismological fittings to a par- 
ticular object. Rather, we concentrate mainly on the con- 
sequences of element diffusion for the chemical profiles ex- 
pected at the DB instability strip. With regard to pulsa- 
tional properties, our attention is directed exclusively to- 
wards the effects of diffusively evolving chemical profiles 
on the adiabatic pulsational properties of the DB white 
dwarf models. In this sense, a detailed description of the 
evolutionary phases prior to the formation of DB white 
dwarfs, particularly regarding the chemical profile for the 
carbon-oxygen core, is not of primary importance for our 
purposes. As for pulsational calculations, we have com- 
puted adiabatic, nonradial g-modes of the DB white dwarf 
models with the pulsational code employed in Corsico et 
al. (2001), appropriately modified to study pulsating DB 
white dwarfs. In particular, the treatment we follow to as- 
sess the Brunt- Vaisala, frequency (N) is that proposed by 
Brassard et al. (1991). 

3. Evolutionary and pulsational results 

We begin by examining Fig. [21 which illustrates the evolu- 
tion of the internal 4 He and 12 C distribution as a function 
of the outer mass fraction q for 0.60-M Q DB white dwarf 
models with a helium content of 8x 10 -3 Af*. In particular, 
the expectation at three selected values of T e s that cover 
the temperature range for pulsating DB white dwarfs is 
depicted. Clearly, diffusion processes appreciably modify 
the chemical profiles as the white dwarf evolves, making 
essentially the bulk of helium float to the surface and heav- 
ier element sink out of surface layers. Note that, as found 
in previous studies (Dehner & Kawaler 1995; Gautschy & 
Althaus 2002 and Fontaine & Brassard 2002) the forma- 
tion of a double-layered structure is expected by the time 
the white dwarf has reached the domain of pulsational in- 
stability. Indeed, diffusion proceeds efficiently, giving rise 
to pure helium outermost layers. Specifically, we find that 
at T cff « 27500K, the O.6-M DB model is characterized 
by a pure helium mantle of 2 x 10 — 6 ikf* (curve b) and 
an underlying intershell (still-uniform) region rich in he- 
lium, carbon and oxygen. We note that the pure helium 
mantles of our models are somewhat larger than those 
reported by Dehner & Kawaler (1995) and Fontaine & 
Brassard (2002). In part, the difference is due to the in- 
clusion of thermal diffusion in our simulations (see sec- 
tion 2). Interestingly, the chemical profile evolves appre- 
ciably with the further evolution through the DB insta- 
bility strip, as shown by Fig. [21 It is worth noting that, 
during such stages, element diffusion not only thickens the 




Fig. 2. Evolution of the abundance by mass of 4 He (solid 
lines) and 12 C (dashed lines) as a function of the outer 
mass fraction (q) for O.6O-M0 DB models having a helium 
content of Mh c = 8 x 10 -3 M*. Starting from a model with 
an initially homogeneous envelope (curves a), following 
models (b, c, d and e) correspond to evolutionary stages 
characterized by T eS = 27500, 23000, 20200 and 10200 K. 
Gravitationally induced diffusion leads to the development 
of a double-layered structure. Note also that the outer 
layer abundance distribution evolves appreciably as the 
white dwarf cools down through the stages where DBs are 
found to pulsate. 



pure helium mantle but also it modifies the shape of the 
chemical profile, a feature which is expected to leave some 
traces in the theoretical period distribution, as compared 
with the case in which the shape of the profile is assumed 
to be fixed during the DB white dwarf evolution across 
the instability strip (see later in this section). 

Although the attention of this paper is focused on the 
evolutionary stage corresponding to variable DBs, we have 
followed the further evolution of some sequences down to 
the domain of the carbon-enriched DQ white dwarfs. The 
resulting chemical profile for the case of a 0.6-Mq DB 
model with a helium content of Mh c = 8 x 10~ 3 M» at 
T eS = 10200K is shown by curve e in Fig. [21 Note that 
even at such advanced stages, diffusion has not engulfed 
the entire helium-carbon-oxygen intershell, and the chem- 
ical profile is still characterized by a double-layered struc- 
ture. In particular, the pure helium mantle extends down 
to a fractional mass depth q rj 5 x 10 _5 M*. Note that, 
despite the total helium content amounting to a fraction 
of « 0.01 of the stellar mass, there is a substantial amount 
of carbon stretching outwards from the intershell region. 
We expect then a significant carbon enrichment in the 
surface layers as a result of convective dredge-up of the 
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Fig. 3. Abundance distribution of 12 C in the outer lay- 
ers at (from bottom to top) T off = 14000, 13300, 11900, 
10700, 10200 and 9200 K. Models correspond to a 0.60- 
Mq DB white dwarf having a helium content of Mn c = 
8 x 10~ 3 Af*. Note that, as a result of convection reaching 
the 12 C diffusive tail, the outer layers become substantially 
12 C-enriched. 



carbon diffusive tail, a process which is generally believed 
to be responsible for the observed carbon in DQs. This 
expectation is indeed borne out by Fig. |3J which shows 
the carbon abundance distribution in terms of the outer 
mass fraction at various T c g values for the O.6O-M0 DB 
white dwarf sequence with a helium mass of 8 x 10~ 3 M*. 
Chemical diffusion causes some carbon from deeper layers 
to diffuse upwards (see Eq. 3) to where convection dredge 
it up to the surface. As cooling proceeds, the base of the 
convection zone moves inwards and eventually reaches a 
maximum depth of 4.2 x 10~ 6 M by T cS « 10000K. The 
depth reached by the convection zone, which at this low 
effective temperature is almost independent of the treat- 
ment of convection, is limited by the low conductive opac- 
ity of the degenerate interior. It is worth noting that the 
presence of carbon (and dredged metals) in the outer lay- 
ers renders the helium plasma less transparent, thus pro- 
viding a not so deep convection zone. Our results sug- 
gest that at T c ff = 10000K, the carbon surface abundance 
reaches 1 log n( 12 C)/n( 4 He) « —4. This abundance far 
exceeds the lowest carbon abundance observed in many 
DQs (see MacDonald et al. 1998). We speculate that DB 



In terms of a single-layered profile, a pure helium man- 
tle of ~ 5 x 10 M* corresponds to a total helium con- 
tent of 10" 3 -10" 4 M«. Thus, our quoted surface carbon abun- 
dance is in rough comparable agreement with the prediction of 
MacDonald et al.(1998) for their single-layered sequences with 
a helium content of about 10 _3 M*. 
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Fig. 4. Same as Fig.[2]but for a helium content of Mh c = 
9 x 10 _4 M*. Starting from a model with an initially homo- 
geneous envelope (curves a), following models (b, c and d) 
correspond to evolutionary stages characterized by T e ff = 
27200, 22900 and 19800 K. Here, the shape of the chemical 
profile is notably modified by diffusion processes. In par- 
ticular, shortly after the red edge of the instability strip 
is reached, the double-layered structure has virtually dis- 
appeared (curves d). 



models less massive than 0.6 Mq evolved directly from 
PG1159 configurations should also exhibit such large car- 
bon surface abundances. In fact, even in this case, settling 
of carbon from the outer layers would take long time, 
and accordingly, a substantial amount of carbon would 
be dredged upwards by convection. Additional evolution- 
ary calculations for such low-mass DB white dwarfs are 
required to place this statement on a solid background. 
Hence, the plausibility of an evolutionary link between 
low-mass PG1159 and DQ stars with low detected carbon 
abundance is not clear if canonical convective dredge-up 
is the source of the observed carbon for such DQs. 

In the case of small helium contents, the shape of the 
double-layered profile is strongly affected by element dif- 
fusion. This is illustrated by Figs. 01 and which show 
the internal He and 12 C distribution as a function of q 
for 0.60-Mq DB white dwarf models with a helium mass 
of 9 x 10~ 4 Af» and 1 x 10 _4 M*, respectively. On the ba- 
sis of these figures, we conclude that the double-layered 
structure is altered by chemical evolution to such a de- 
gree that a single-layered configuration emerges even dur- 
ing the stages of pulsation instability. This result is not 
surprising in view of Eq. (3). Indeed, for small helium 
contents, electron degeneracy occurs below the helium- 
carbon-oxygen mantle, and as a result, the second term in 
Eq (3) is dominant, leading to the formation of a single- 
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Fig. 5. Same as Fig.|21but for a helium content of Mhc = 
1 x 10~ 4 M*. Starting from a model with an initially homo- 
geneous envelope (curves a), following models (b, c and d) 
correspond to evolutionary stages characterized by T e s = 
27300, 23600 and 20600 K. Here, the chemical profile bears 
no clear signature of a double-layered structure even at the 
hottest T c ff values considered (curves b). Models, instead, 
are characterized by a single-layered structure. 



chemical structure in DB white dwarfs 




Log(q) 

Fig. 6. Same as Fig. 0] but for a stellar mass of 0.75 
Mq . Starting from a model with an initially homogeneous 
envelope (curves a), following models (b, c and d) cor- 
respond to evolutionary stages characterized by T e g = 
27200, 23100 and 19900 K. Note that, a single-layered pro- 
file emerges at larger T e ff values as compared with the 0.6 
Mq model. 



layered structure in relatively short ages. Specifically, for 
Mjie = 9 x 10 _4 M* a single-layered configuration occurs 
by the time the star reaches the end of the instability 
strip, and for Mh c = 1 x 10~ 4 M„, a single-layered con- 
figuration sets in even by the onset of pulsation instabil- 
ity. Hence, if post-AGB, DB white dwarf progenitors are 
formed with a helium content smaller than w 10 -3 M*, 
the double-layered structure is expected to become single- 
layered by the time the star reaches the red edge of the 
DB instability strip. 

The dependence of the chemical profile evolution upon 
the stellar mass is shown in Fig. El to [HI In particular, in 
Fig. H3 we compare the evolution of the chemical profile 
for 0.75-M© DB models with the same helium content 
than that considered in Fig.H(M He = 9 x 10" 4 Af»). The 
role of diffusion in relatively massive DB white dwarfs is 
clearly documented in this figure. Indeed, here, over most 
of the temperature range of DB variables, the chemical 
profile bears no signature of a double-layered structure. By 
contrast, in the case of 0.6-M© models the double-layered 
profile turns into a single-layered one only after evolution 
has proceeded to effective temperatures somewhat below 
the red edge of the instability strip. In Figs. and |5] we 
compare the shape of the helium profile for various stellar 
masses at a fixed T e g . Note that the shape of chemical 
profile is very dependent on the stellar mass of the star. 



Next, we examine the pulsational properties of our 
models. We shall consider non-radial (I — 1) g-modes with 
periods Pk (k being the radial overtone) in the range of 
100-1000 s. Figuresl51ll0landlllldisplav the forward period 
spacing (APk = Pk+i — Pk) as a function of period (pe- 
riod spacing diagrams) for the case of 0.6-M Q DB white 
dwarf models with helium contents of Mh c = 8 x 10 _3 Mh,, 
9 x 10~ 4 M* and 1 x 10 _4 M*, respectively. In each figure 
results corresponding to a sequence of four selected mod- 
els at decreasing effective temperatures (as labeled in left 
panels) covering the pulsational strip are depicted. The 
left column of panels (solid lines) depicts the APk dis- 
tributions as predicted by calculations in which element 
diffusion has been accounted for, and the central column 
of panels (dashed lines) shows the same information but 
for the case in which diffusion processes have been locked 
in the evolutionary computations, e.g., the chemical pro- 
files remain fixed with further evolution from a T Q g value 
corresponding approximately to the blue edge of the DB 
instability strip. Finally, the right column of panels docu- 
ments the differences between the values of APk as given 
by diffusively evolving profiles and the fixed ones. 

An inspection of the figures indicates that in both sets 
of computations the pulsation mode spectrum changes 
markedly as the white dwarf models cool down to the 
red extreme of the DB instability strip. In the case of 
models with diffusion, such changes are caused by: (1) the 
changes the mechanical and thermal structure of models 
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Fig. 7. Abundance by mass of 4 He as a function of the 
outer mass fraction (q) for 0.60, 0.75 and 0.85-M Q DB 
models (solid, dashed and dotted lines, respectively), hav- 
ing a helium content of Mjj_ e = 8 x 10~ 3 M*. Models cor- 
respond to T c ff = 23900K. 
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Fig. 8. Same as Fig.Qbut for a helium content of Mh c = 
9 x 10 _4 M*. Solid, dashed and dotted lines correspond to 
DB models with 0.60, 0.70 and O.75-M at T cff = 23700K. 

in response to energy looses, changes which are translated 
mainly into a decreasing in the Brunt- Vaisala frequency 
at the core region and the consequent increasing of pe- 
riods and period spacings, and (2) the evolution of the 
chemical profiles caused by diffusive processes, thus giving 
rise to changes in the location of the well known Brunt- 



Vaisala frequency's bumps associated with chemical inter- 
faces. Instead, in the case of models in which chemical dif- 
fusion is not considered, the changes in the period pattern 
are due only to the secular changes resulting from cool- 
ing (item 1), e.g., the Brunt- Vaisala frequency decreases 
at the innermost regions but the profile of the Ledoux 
term B (a quantity closely involved in the computation 
of the Brunt- Vaisala frequency at the chemical interfaces; 
see Brassard et al. 1991 for details) remains almost fixed. 
Interestingly enough, during the instability strip the pe- 
riod spacing distribution is strongly modified according to 
whether element diffusion is considered or not. 

It is interesting to note that the period spacing dia- 
grams of the pulsation mode spectrum exhibits a clear 
mode-trapping substructure when a multi-layered struc- 
ture is present, and more importantly how this substruc- 
ture change due to diffusion. In fact, from left panels of 
Fig. it can be seen that there is substructure at all evo- 
lutionary stages depicted, whereas in Fig. 1101 the period 
spacing diagram becomes steadly simpler as the T e g de- 
creases. Finally note that in Fig. El the period spacing 
diagram is characterized by a fairly clean structure at all 
temperatures shown. The occurrence of substructures in 
the period diagrams is of course intimately related to the 
presence or not of the double-layered structure, a fact that 
can be related directly to Figs. \M and El That is, the 
substructure is determined by whether or not the pure 
helium layer has had time to engulf the whole helium- 
carbon-oxygen intershell. In the case of a single-layered 
structure there is a single strong mode-trapping cycle, oth- 
erwise there are two cycles of differing strength. 

Finally, the period spacing differences between model 
with and without diffusion (as they were defined before) 
are by no means negligible and, as can be visualized from 
the plots, amounts to « ±8 s for all the helium contents we 
have considered. These differences far exceed the quoted 
values of rms differences between the observed and cal- 
culated period spacings (cap) in current detailed aster- 
oseismological fits to pulsating DB white dwarf stars. In 
light of these results, we conclude that the dependence of 
the shape of the evolving chemical profiles on the effective 
temperature cannot be entirely ignored, that is, in order 
to maintain physical consistency, these parameters cannot 
be considered as independent ones when attempts at per- 
forming detailed seismological fit to a particular star are 
made. These parameters are actually not coupled if the 
helium mass is less massive enough for a single-layered 
profile to emerge. 

4. Conclusions and discussion 

In this paper we have explored the evolution of white 
dwarf stars with helium-rich atmospheres (DB) in a self- 
consistent way with the predictions of time-dependent el- 
ement diffusion. The focus is on the chemistry variations 
along the effective temperature range where these stars 
are found to pulsate. There exist theoretical and observa- 
tional evidence (Dehner & Kawaler 1995; Dreizler & Heber 
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Fig. 9. Period spacing diagrams (for I = 1 modes) for 
the case of 0.6-M Q WD models with a helium content of 
Mjie = 8x 10 _3 M*. From top to bottom panels correspond 
to decreasing effective temperatures (labels in the left pan- 
els). The left (centre) column of panels shows the results 
obtained when element diffusion has (not) been considered 
in the computations. Right panels depict the differences 
between the period spacing as given by diffusively evolv- 
ing profiles and a fixed one predicted by diffusion at T e g w 
27000K. 

1998) suggesting that the majority of the DB white dwarfs 
could be the descendants of PG1159 stars, the hot and 
hydrogen-deficient remnants of post-AGB stars that ex- 
perienced a very late helium thermal pulse on their early 
cooling phase (see Herwig et al. 1999). In this regard, we 
have assumed for our starting DB white dwarf models 
an uniform envelope chemical composition rich in helium, 
carbon and oxygen. For this study, we have considered 
DB white dwarf models with stellar masses of 0.60, 0.65, 
0.70, 0.75 and 0.85 M . The mass of the helium enve- 
lope (Mho) was varied in the range 10~ 2 — 10 _4 M*. For a 
consistent treatment with the predictions of element dif- 
fusion, OPAL radiative opacities for arbitrary metallicity 
including carbon-and oxygen-rich compositions have been 
employed. Our treatment of element diffusion, based on 
the multicomponent gas description of Burgers (1969), in- 
cludes the effects of thermal diffusion, which, as we have 
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Fig. 10. Same as Fig.lSJbut for a helium content of Mh 6 = 
9 x 10- 4 M*. 

shown, lead to DB models with thicker pure helium man- 
tles. 

In view of the lack of evolutionary calculations with 
time-dependent diffusion for DB white dwarfs with differ- 
ent stellar masses and helium envelopes, an initial ques- 
tion to be addressed in this work was whether the double- 
layered chemical profile expected for hot variable DBs (as- 
suming that these stars descend from PG1159 stars) could 
be markedly modified by further chemical evolution be- 
fore the red edge of the instability strip is reached. This 
is an important point because as recently demonstrated 
by Fontaine & Brassard (2002) the presence of a double- 
layered structure leads to a distinct theoretical period dis- 
tribution from that predicted by single-layered configura- 
tions. Within the scope of the present paper is also the 
study of the consequences of diffusively evolving chemi- 
cal profiles for the adiabatic pulsational properties of DB 
white dwarfs. 

Our results show that, as found in previous studies, a 
double-layered envelope is expected by the time the white 
dwarf reaches the domain of the DB pulsational instability. 
In particular, element diffusion not only thickens the pure 
helium mantle as cooling proceeds but also it modifies the 
shape of the chemical profile, an aspect to be considered in 
asteroseismological fittings of pulsating DB white dwarfs. 
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Fig. 11. Same as Fig.Elbut for a helium content of Mr c = 
1 x l(r 4 M„. 



In this regard, we find that the diffusively evolving shape 
of the envelope composition translates into a distinct be- 
haviour of the theoretical period distribution as compared 
with the case in which the shape of the profile is assumed 
to be fixed during the DB white dwarf evolution across the 
instability strip. In particular, we find that the presence 
of a double-layered structure causes the period spacing- 
diagrams to exhibit mode-trapping substructures which 
change due to diffusion. 

Another finding of this work is the fact that the double- 
layered structure is altered by further chemical evolution 
to such a degree that, depending on the helium content, 
a single-layered configuration emerges even during the 
stages of pulsation instability. In particular, for a 0.6 Af© 
DB white dwarf with Mn e = 10~ 3 Af* this occurs by the 
time the star reaches the T e g value of « 20000K, and for 
Mne = 10 _4 M*, a single-layered configuration is expected 
at T eff ss 27000K. Thus, if post-AGB, DB white dwarf pro- 
genitors are formed with a helium content smaller than 
w 10~ 3 A/*, the double-layered structure is expected to 
become single-layered by the time the star reaches the red 
edge of the DB instability strip. A helium envelope mass of 
« 10 _3 M* is more appropriate for a massive white dwarf. 
In particular, we find that the double-layered profile in a 



w 0.8 Mq DB turns into a single-layered one well before 
the red edge of the DB instability domain is reached. 

The existence of 0.6 Mq DB white dwarfs with a he- 
lium content less massive than, say, « 10 _3 Af* is difficult 
to reconcile with the predictions of the born-again scenario 
for the formation of hydrogen-deficient post-AGB stars. 
Indeed, Herwig et al. (1999) find that after the end of the 
late helium thermal pulse the total helium content left in 
the star amounts to A/hc ~ 10 _2 Af». Now, there exist ob- 
servational evidence suggesting that post-AGB mass-loss 
episodes could reduce (and even remove) the mass of the 
helium-rich envelope considerably. As recently emphasized 
by Werner (2001) it is possible that the surface chemistry 
of all PG1159 stars is determined not only by previous 
interior mixing processes but also by mass-loss episodes. 
Mass-loss rates in the range 10~ 7 — lO _8 M0/yr are de- 
tectable in many luminous PG1159 stars. In addition, ten- 
tative evidence for the persistence of mass-loss rates of 
the order 10~ 7 — 10 -10 M Q /yr down to the domain of hot 
helium-rich white dwarfs has been presented (see Werner 
2001). The existence of such mass-loss rates would imply 
that most of the helium-rich envelope of DB progenitors 
could be substantially reduced during the 10 5 — 10 6 yr 
elapsed during the post-AGB stage. It is worth mentioning 
that the existence of PG1159 stars with a helium content 
as low as 10 -3 Mq has been suggested by asteroseismol- 
ogy in at least one of these stars with a stellar mass of 0.6 
Af Q (Kawaler & Bradley 1994), thus implying the occur- 
rence of modest mass-loss during evolution to the PG1159 
phase. In light of such evidence, the possibility of the exis- 
tence of some DB white dwarfs with helium contents lower 
than 10~ 3 A/* could not be discounted. According to our 
calculations, a double-layered structure would thus not be 
expected for such white dwarfs once they have reached the 
domain of DB pulsational instability. 

Finally, we have extended the scope of our evolutionary 
calculations down to effective temperatures characteristics 
of the helium-rich carbon-contaminated DQ white dwarfs, 
the supposed cooler descendants of DBs. The presence of 
carbon in the surface of such stars is widely believed to 
result from convective dredge-up of the carbon diffusive 
tail by the superficial helium convection zone (Pelletier 
et al. 1986, see also Koester et al. 1982). New dredge- 
up carbon simulations in DQ stars by MacDonald et al. 
(1998) suggest that the low carbon abundance observed 
in many DQs could be explained in terms of a dredge- 
up effect only if such white dwarfs are characterized by 
very thick helium envelopes. As our results shown, if we 
assume that DQ white dwarfs may be linked to the post- 
AGB PG1159 stars via the PG1159-DB-DQ connection, 
then a 0.6 Mq DQ white dwarf with a helium content 
more massive than 10 _2 Af» would be characterized by a 
double-layered structure. In particular, we find that the 
pure helium mantle extends down only to a fractional 
mass depth g«5x 10~ 5 Af» at T e $ = 10000K, shallow 
enough for carbon to be dredged-up to the surface with 
abundances far exceeding the observational predictions. 
If canonical convective dredge-up theory is correct, our 
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results seem to discount an evolutionary connection be- 
tween some DQs, those having low carbon surface abun- 
dances, and the PG1159 stars. Full evolutionary calcula- 
tions taking into account the evolutionary stages leading 
to the white dwarf formation, and particularly the chem- 
ical structure emerging from the born-again scenario, are 
required to place these assertions on a more solid basis. 
Work in this direction is in progress. 

Detailed tabulations of the chemical 
profiles of our models are available at 
htt p : / / www . fcaglp . unlp . edu . ar / evolgroup / 
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